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Abstract. Let X be the branching particle diffusion corresponding to the 
operator Lu + /3(u 2 — u) on D C R d (where j3 > and 0^0). Let A c denote 
the generalized principal eigenvalue for the operator L + (3 on D and assume 
that it is finite. When A c > and L+/3— A c satisfies certain spectral theoretical 
conditions, we prove that the random measure exp{— X c t}Xt converges almost 
surely in the vague topology as f tends to infinity. This result is motivated by 
a cluster of articles due to Asmussen and Hering dating from the mid-seventies 
as well as the more recent work concerning analogous results for superdiffusions 
of [ETllEW] , We extend significantly the results in |AH76llAH77l and include 
some key examples of the branching process literature. As far as the proofs 
are concerned, we appeal to modern techniques concerning martingales and 
'spine' decompositions or 'immortal particle pictures'. 
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1. Introduction and statement of results 

1.1. Model. Write C Z ' V (D) to denote the space of i times (i = 1,2) continuously 
differentiable functions with all their ith order derivatives belonging to C n (D). 
[Here C T? (D) denotes the usual Holder space.] Let D C M. d be a domain and 
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consider Y = {Y t ; t > 0}, the diffusion process with probabilities {P^, x G D} 
corresponding to the operator 



where the coefficients aij and hi belong to C 1 ' 7 * , i,j = 1 , d, for some 77 in (0, 1], 
and the symmetric matrix a = {ai.j} is positive definite for all x G D. At this point, 
we do not assume that Y is conservative, i.e. Y may get killed at the Euclidean 
boundary of D or run out to infinity in finite time. 

Furthermore let us first assume that < (3 G C V (D) is bounded from above on D 
and (3 ^ 0. The (strictly dyadic) (L, (3; D)-branching diffusion is the Markov process 
with motion component Y and with spatially dependent rate (3, replacing particles 
by precisely two offspring when branching and starting from a finite configuration of 
individuals. At each time t > 0, the state of the process is denoted by X t G M. (D) 
where 



We will also use the following notation: X = {X t : t > 0} has probabilities {P^ : \i G 
M (D)}, and is expectation with respect to P^. As usual, (/, fi) := J D f{x) (i(dx) 
and (/, g) := f D f(x)g(x)dx, where dx is Lebesgue measure, and so (f,gdx) = 



When (3 is not bounded from above, one may wonder if the (L, (3; Z))-branching 
diffusion is still well defined, in particular, whether the global (or even local) mass 
may blow up in finite time. To treat the case with /3's which are not upper bounded 
we will need to consider more general branching diffusions. For a 1 weighted branch- 
ing diffusion', the particles do not necessarily carry unit mass. At each time t > 0, 
the state of the process is X t e M (D) where 



Next we need a definition. Let 

A c = X C (L + [3, D) := inf{A G M : 3u > satisfying (L + (3 - X)u = in D} 

denote the generalized principal eigenvalue for L + (3 on D. By standard theory, 
A c < 00 whenever (3 is upper bounded and, for general (3, there exists an h > 
satisfying that (L + (3 — X)h = 0, whenever A c < 00. 

In the latter case, let us define the M (D)-valued process W as follows. Each 
particle performs a Y motion and carries weight h(x) at x G D, and furthermore, 
when the particle's clock rings, according to the spatially varying rate (3, the par- 
ticle splits into two offspring, which perform independent Y motions and carry 
weights according to the function h, etc. Since (3 is not upper bounded, we have 
to rule out finite time explosions. Fortunately, since h is a 'harmonic' function, 
it is a straightforward exercise to show that the total mass process \W\ is a su- 
permartingale, and in particular, it is a.s. finite for all t > 0. Then, since X is 



well defined, so is the M. (D)-valued process X defined by := e Xat h 1 , (i.e. 
X t {B) := e^ih-i-lB,^), t>0, Be D Borel). 



Therefore, from now on, we relax the assumption that swp D {3 < 00 and replace 
it with the much milder assumption A c < 00. 



(1) 



L = -V-aV + 6-V on R d 




(fg,dx) = (f,g) D 
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1.2. Motivation. This paper concerns growth of mass on compact domains of 
branching particle diffusions. In doing so we address a gap in the literature dating 
back to the mid-seventies when the study of growth of typed branching processes 
on compact domains of the type space was popularized by Asmussen and Hering. 
Also we complement a recent revival in this field which has appeared amongst the 
superprocess community. 

Before discussing main results, we shall introduce the topic in detail. 

Definition 1 (Local extinction). Fix {i E M (D). We say that X exhibits local 
extinction under P^ if for every Borel set B CC D , there exists a random time Tg 
such that 

P^Ijb < oo) = 1 and P^(X t (B) = for all t > t b ) = 1. 

[Here B CC D means that B is bounded and its closure is a subset of D.] 

Local extinction has been studied by |P96j . [EPj (for superprocesses) and |EK 
(for branching diffusions). To explain their results, recall that we assume that the 
generalized principal eigenvalue for L + (3 on D is finite. In fact, A c < if and 
only if there exists a function h > satisfying (L + (3)h = on D - see Section 4.4 
in |P95j . Following the papers P96, EP] where similar issues were addressed for 
superprocesses, in |EK] the following was shown. 

Theorem 2 (Local extinction versus local exponential growth). Let ^ /i e 
M(D). 

(i) X under P^ exhibits local extinction if and only if there exists a function 
h > satisfying (L + (3)h = on D, that is, if and only if A c < 0. 

(ii) When A c > 0, for any A < A c and ^ B CC D, 

P M (limsup tToo e- xt X t {B) = oo) > 0, 

iyiimsup nco e~ Xct X t (B) < oo) = 1. 

In particular, local extinction/local exponential growth does not depend on the initial 
measure fi G A4 (D) . 

(In [EKj it is assumed that (3 is upper bounded, whereas in [EPj only the finite- 
ness of A c is assumed. The proofs of [EKj go through for this latter case too.) 
On closer inspection this last theorem says that when A c < mass 'escapes out 
of £?' even though the entire process may survive with positive probability. (If 
Y is conservative in D for example then it survives with probability one) . Further, 
when A c > mass accumulates on all nonempty bounded open domains and in 
such a way that with positive probability this accumulation grows faster than any 
exponential rate A < A c . On the other hand, mass will not grow faster than at 
the exponential rate A c . It is natural then to ask whether in fact A c gives an exact 
growth rate or not. That is to say, for each ^ B CC D do the random measures 
{exp{— A c i}A t : t > 0} converge in the vague topology almost surely? Further, 
can one identify the limit? This is precisely the object of interest of a variety of 
previous studies for both branching diffusions and superprocesses which we shall 
now review. 

We note already here that the process in expectation is given by the linear 
kernel corresponding to the operator L + (3 on D. Therefore, trusting in SLLN for 
branching processes, one should expect that the process itself grows like the linear 
kernel too. On the other hand, it is easy to see that the linear kernel does not 
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in general scale precisely with exp{— X c t} but rather with f(t) cxp{— A c <}, where 
/ grows to infinity as t — > oo and at the same time is subexponential. (Take for 
example L = A/2 and /3 > constant on R d , then f(t) = t d/2 .) In fact the 
growth is pure exponential if and only if L + (3 is product-critical (see later in this 
subsection). Proving SLLN seems to be significantly harder in the general case 
involving the subexponential term /. 

In the late seventies Asmussen and Hering wrote a series of papers concerning 
weak and strong laws of large numbers for a reasonably general class of branching 
processes which included branching diffusions. See AH76] and [AH77] . In the 
context of the branching diffusions we consider here one can summarize briefly their 
achievements by saying that, when D is bounded, for a special class of operators 
L + (3, the rescaled process {exp{— X c t}X t : t > 0} converges in the vague topology, 
almost surely for branching diffusions. Further, for the same class of L + /3 when D 
is unbounded they proved that there exists the limit in probability of exp{— X c t}X t 
as t | oo (in the vague topology). The class of L+(3 alluded to they called 'positively 
regular'. The latter is a subclass of the class V*(D) (the class that we shall work 
with) given below. 

A more detailed comparison with |AH76[ IAH77] as well as the discussion on 
related results on superprocesses is deferred to Section [2J 

Before we give the definition of the basic classes of operators that we shall use, 
V P (D) and V*(D), we need to recall certain concepts of the so-called criticality 
theory of second order operators. The operator L + (3 — \ c is called critical if the 
associated space of positive harmonic functions is nonempty but the operator does 
not possess a (minimal positive) Green's function. In this case the space of positive 
harmonic functions is in fact one-dimensional. Moreover, the space of positive 
harmonic functions of the adjoint of L + /3 — A c is also one dimensional. 

Assumption 3. Suppose we choose representatives of these two spaces to be <fi 
and <p respectively. Throughout the paper and without further reference, we will 
always assume that L + (3 — A c product-critical, and in this case we pick 4> and 4> 
with the normalization (4>,4>) = 1. 

We now define the classes V p (D) and V* (D) . Since we want to talk about 
spatial spread on a generic domain D, we fix, for the rest of the paper, an arbitrary 
family of domains {D t , t > 0} with D t CC D, D t f D. (For D = R d , D t can be 
the i-ball, but we can take any other family with D t CC D, D t /\ D too.) 

Definition 4 (V p (D) and V* (£>)). For p > 1, we write L + (3 £ V v (£>) if 

(i) : A C = A C (L + /?;£>) > 0, 

(ii) : {<j} p , (f>) < oo, in which case we say that L + j3 — A c is product p-critical. 

Let q(x, y, t) be transition density of L+(3 and Q(x, y, t) := q{x, y, t)—e Xat (f)(y)(j){x). 
We write L + (3 £ V* (D) when the following additional conditions holds. 

(iii) : For all 7^ B CC D there exists a positive function £ with £(t) f 00 as 
1 1 00 and such that 

\Q{x, y, C(t))| < ot t 4>{y)ct>{x), x£D u y£B, t>0, 

where limtjoo e~ Xct at — 0. 



STRONG LAW OF LARGE NUMBERS FOR BRANCHING DIFFUSIONS 



5 



(iv): Furthermore, for any given x G D there exists a function a : [0, oo) — > 
[0, oo) such that ((at) = 0(t) as t — > oo and for all 6 > 0, 

(3n ,Vn > n : supp(X„ 5 ) C D anS ) = 1. 

Let p(x,y,t) denote the transition density of the diffusion corresponding to the 
operator (L + (3 — A c )^. Then p(x, y, t) — e _Act </>(j/)0 _1 (x)g(a;, y, i), and thus, (hi) 
is equivalent to 

(iii*): With the same C, as in (hi), 

p(x,y,C(t)) _ 1 =Q 
#(y) 

Some heuristics to help find suitable a and C will be discussed in Section [3] 

Remark 5 (Ergodicity). Note that criticality is invariant under ^-transforms. 
Moreover, an easy computation shows that <p and cj> transforms into 1 and <jxp 
respectively when turning from (L + (3 — A c ) to the /i-transformed (h — <fi) opera- 
tor (L + P — X c )^ = L + acj) -1 V0 • V. Therefore product criticality is invariant under 
/i-transforms too (this is not the case with product p-criticality when p > 1). Fur- 
ther, for operators with no zeroth order term, it is equivalent to positive recurrence 
(ergodicity) of the corresponding diffusion process. In particular, (L + /3 — A c )^ cor- 
responds to an ergodic diffusion process provided (L + ft — A c ) is product critical 
(see [P95] . Section 4.9). o 

1.3. Main results. With the following theorem we wish to address the issue of 
almost sure convergence in the vague topology of {exp{— X c t}X t : t > 0} for branch- 
ing diffusions with L + j3 belonging to V*(D), p > 1 thus generalizing the results 
of Asmussen and Hering. 

Note that since L + f3 — A c is critical, <fi is the unique (up to constant multiples) 
invariant positive function for the linear semigroup corresponding to L + (3 — X c 
(Theorem 4.8.6. in |P95j ). Let {St}t>o denote the semigroup corresponding to 
L + (3. It is a standard fact (sometimes called 'the one particle picture') that 

(2) S t (g)(x) = E 5x (g,X t ) 

for all nonnegative bounded measurable g's. Even though (f> is not necessarily 
bounded from above, St(</>) makes sense and @ remains valid when g is replaced 
by <p, because <fr can be approximated with a monotone increasing sequence of g's 
and the finiteness of the limit is guaranteed precisely by the invariance property 
of (f>. By the invariance of 0, Es D .e~ Xct ((f), X t ) — e _A °* St{<j>)(x) = <p(x), which is 

sufficient together with the branching property to deduce that = {wf; t > j 

is a martingale where 

Wf = e- x ° t {4>,X t ), t>0. 

Indeed note that 

E Sm (e-^+^^Xt+s)]^ = e^E Xt (e- x < s (^X s )) = e' x ^{^X t ). 

Being a positive martingale, Ps^ -almost sure convergence is guaranteed, and the 
a.s. martingale limit := lim^oo Wf appears in the following main conclusion. 



lim sup 
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Theorem 6 (SLLN). Assume that L + [3 e V*{D) for some p E (1,2] and 
(I3(j) v ,4>) < oo. Then, 

(3) lim tToo e-^ig, X t ) = (g, j)W+, g G C+ (£>) 

ftoZds Pa x - a.s. /or a// x e D, and Eg x (W*,) = 1. 

Moreover, if sup D (3 < oo then the restriction p € (1,2] can 6e replaced by p > 1. 

We close this subsection with the Weak Law of Large Numbers. Here we change 
the class V*{D) to the larger class V P (D) and get convergence in probability instead 
of a.s. convergence. It is important to point out, however, that the class V*(D) is 
already quite large — see Section [3j where we verify that key examples from the 
literature are in fact in V* (D) and thus obey the SLLN. 

Theorem 7 (WLLN). Suppose that L + (3 e V P {D) for some p £ (1,2] and 
(f3(j> p , <j>) < oo. Then for all x € D, holds in Ps x - probability and Es x (W^) = 1. 
Moreover, if sup D (3 < oo then the restriction p 6 (1,2] can be replaced by p > 1. 

1.4. Outline. The rest of this paper is organized as follows. In Section [2] we embed 
our results into the literature, while in Section [3] we discuss some key examples for 
the SLLN. The proofs are given in Section 4. 

2. Detailed comparison with some older results 

The methods of Asmussen and Hering were based for the most part on classical 
techniques of truncation and applications of the Borel-Cantelli Lemma. Using this 
method, they proved the convergence of e~ Xct (X t , g) for all < g € L 1 (cf)(x)dx) . 
It is also worth noting that the generic strength of their method extended to many 
other types of branching processes; discrete time, discrete space and so on. 

Interestingly, preceding all work of Asmussen and Hering is the single article fW] 
(later improved upon by [B]). Watanabe demonstrates that when a suitable Fourier 
analysis is available with respect to the operator L+(3, then by spectrally expanding 
any g € (D) , the space of nonnegative, continuous and compactly supported 
functions, one can show that {(g,X t ) : t > 0} is almost surely asymptotically 
equivalent to its mean. From this the classic Strong Law of Large Numbers for 
dyadic branching Brownian motion in M. d is recovered. Namely that when L = A/2 
and [3 > is a constant, 

lim tToo t d / 2 e- 0t X t (B) = {2irf 2 \B\ x 

where B is any Borel set, \B\ is its Lebesgue measure and is a strictly positive 
random variable depending on the initial configuration /i 6 Ai(R d ). The operator 
1/2A + /3 does not fall into the class Vi(D). For an analogous result on supercritical 
super-Brownian motion see [Ej . 

Let us discuss now how our assumptions relate to the assumptions imposed in 
the article |AH76j . 

In [AH76] the domain is bounded (and even one dimensional) when the Strong 
Law of Large Numbers is stated for branching diffusions; on general domains, only 
convergence in probability was obtained. Furthermore, in [AH76 the notion of 
positively regular operators was introduced. In our context it first means that 

(A) A c > (they call this property 'supercriticality'). 

(B) is bounded from above, 
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(C) (<M> <oo. 



Obviously, (B-C) is stronger than the assumption (0, (f>) < oo (product-criticality) . 

Secondly, {S t }t>o, the semigroup corresponding to L + (3 (the so called 'expec- 
tation semigroup') satisfies the following condition. If r/ is a nonnegative, bounded 
measurable function on M. d , then 

(D) S t (r])(x) = (rj, </>) <j>(x) [e Xct + o (e Act )] as t f oo, uniformly in r). 

Let T t be the semigroup defined by T t (f) := Sf{f) = \S t {4>f), for all < / 
measurable with (f>f being bounded. Then T t correspond to the /i-transformed 
(h = <j>) operator Lq. Recall that Lq corresponds to a positive recurrent diffusion 
process. Then, assuming that 4> IS bounded, it is easy to check that the following 
condition would suffice for (D) to hold: 



where || • |) denotes sup-norm. However this is not true in most cases on unbounded 
domains (or even on bounded domains with general unbounded coefficients) because 
of the requirement on the uniformity in x. (See our examples in Section[3] — neither 
of the examples on M d satisfy (HI).) 

Turning to superprocesses, there would seem to be considerably fewer results 
of this kind in the literature (see the references [Dl |Dy| [E] for superprocesses in 
general). The most recent and general work in this area we are aware of are |ET| 



In |ET] it was proved that (in the vague topology) {exp{— X c t}X t : t > 0} 
converges in law where X is the so called (L, (3, a, R d )-superprocess (with a being 
the 'activity parameter') satisfying that L + (3 S V (D) and that acj) is bounded 
from above. (An additional requirement was that (cf>, /i) < oo where /i = X is the 
deterministic starting measure. ) The long and technical proof relied heavily on the 
theory of dynamical systems applied to the Laplace transforms of {e~ Xct cf)Xt, t > 0}. 

In EW the convergence in law was replaced by convergence in probability. 
Furthermore, instead of M. d a general Euclidean domain D C M. d was considered. 
The heavy analytic method of |ETj was replaced by a different, simpler and more 
probabilistic one. The main tool was the introduction of a 'weighted superprocess' 
obtained by a 'space-time 77-transform'. 



In this section we give examples which satisfy all the assumptions we have, and 
thus, according to Theorem [51 obey the SLLN. (Those examples do not fall into 
the setting in [AH76llA"H77] .) 

Before we turn to the specific examples, we give some heuristics. Although these 
are not actually needed for understanding the examples, we feel that the reader 
'gets a more complete picture' by first reading them. 

Remark 8 (Expectation calculations and local vs. global growth rates). From {5}, 
we have 



(4) 



*T°° xt zD \\g\\<l 



limsup sup (g, <j>4>) 1 T t (g) - (g , <j>4>) = 0, 



EW, E . 



3. Examples 



E x (l{ ■ ed y },X t ) — e 



.A e t 



<Kv) 



s 
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and then, by ergodicity, 

e-^E x {l,X t ) = <t>{x) [ P{t ^ ] dy -> 4>{x) [ 4>{y) dy, as t -> oo. 

Hence, if (cf>, 1) < oo, then the global population growth is the same as the local 
population growth, whereas, if (cf>, 1) < oo the global growth rate exceeds the local 
growth rate. 

Remark 9 (Heuristics for a t and C(i)). One may wonder how one can find a 
function a as in Definition [5Jiv) . In fact, this will often be straightforward to find. 

Suppose that Y is conservative and fix x £ D. (If Y is not conservative, then 
there is no function a satisfying (iv). Indeed, D t f implies that U«r-Dt C Dt, 
for T > 0. Now, if Yr = A with positive probability for some T > 0, then the 
requirement D aT CC D cannot hold.) If, for example, we can pick a deterministic 
increasing function a such that, for all S > 0, 

oo 

J~] Ps x (supp(X n5 ) £ D a „J < oo, 

n=l 

then Borel-Cantelli says that the function a is an appropriately choice. Since the 
probability one particle is present in a set is trivially dominated by the expected 
numbers in that set, it will be much easier to check that 

oo 
n=l 

If we can choose at such that, for some e > 0, 

/ %^ <e -(Ac + e)t 

y| v i>a t 4>{y) 

then we will have satisfied 

oo 

E x {l D c ,X nS ) < oo. 



OO 

E 



Heuristically, if the spine transition density p(t, x, y) converges to its equilibrium 
4>{y)4>{y) sufficiently quickly even for very large y, we might hope to take 

a t « 4>- 1 (e-^ 4 ) 

If the spine starts at a very large position, since it is ergodic it will tend to move 
back toward the origin and Ventcel-Friedlin large deviation theory suggests that it 
will 'closely' follow the path of a deterministic particle with the same drift function. 
We can use this to guess for a suitable form for At least heuristically, to find 
out how far away the spine particle may start in order that it both returns to the 
vicinity of the origin and then ergodizes towards its invariant measure before large 
time t, we can solve the deterministic differential equation 

when L = ha(x)A — fi{x) ■ V, and take a little larger than \f{t)\. 

Indeed, these heuristics appear to the correct form for both at and £(t) in the 
examples considered below. 
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Example 10 (OU process with quadratic breeding rate). Let a,fi,a,b > and 
consider 

L := -a 2 A - fix ■ V on R d 

corresponding to an (inward) Ornstein-Uhlenbeck process and let (3(x) := bx 2 + a. 
Since L corresponds to a recurrent diffusion and is a smooth function with /3 > 
and P ^ 0, it follows that A c > (see Chapter 4 in [P95J). The equilibrium 
distribution for L is given by a normal distribution, 

^ )= (i) d/2eXP {-|^ 2 }- 
Suppose that /i > a\f 7 2h. Defining 7 ± :— {^i ± \J \i 2 — 2ba 2 ^j , for the principle 
eigenvalue problem with (L + /?)</> = \ c <fr we can take 

A c := <t 2 7~ + a, <j)(x) :— c~ exp{7~x 2 } and (f>(x) — c + exp{— "f + x 2 } 

where c" := (l - {2bcr 2 / fi 2 ))^ , c+ := c~ (^/{tkj 2 )) * . Note that 4>{x) := 4>(x)n{x) 
and L + [5 is a self-adjoint operator with respect to tt. 

Some calculations using the 'one-particle picture' (equation [2]) reveals that, in 
expectation, the support of the process grows like \/X c t/j+ and one can pick a(t) = 
yAt/7+ for any A > A c and condition (iv) in Definition [4] will hold. 

The spine is also an (inward) Ornstein-Uhlenbeck process with parameter a :— 
[i - 2-f-<7 2 = vV 2 - Zba 2 with 

(L + (3 - A c )* = L + cr 2 — • V = —a 2 A - ax ■ V on R d , 
cf> 2 

and transition density 

»Eti(^"^' (a/a2)t ) 2 " 

_ e -2(a/<r 2 )t) 

We see that the drift of the inward OU reduces the influence of any starting position 
exponentially in time. Indeed, one can take Ci x ) = (1 + e )( fj2 / Q logs for any e > 
for condition (in*) in Definition [4] to hold. Finally, we trivially note that £(a(t)) = 
0(t) (in fact, only \ogt growth), hence all necessary conditions are satisfied for our 
strong law theorem to hold. 

Note, a strong law for a generalization of this model can be found in [H] where 
the convergence is proved using a martingale expansion for continuous functions 
g € L 2 (tt) (rather than compactly supported g). Almost sure asymptotic growth 
rates (and a.s. support) for the same model are studied in |GHH| . 

This is certainly a non-trivial model and it highlights the strength of our general 
result. In particular, a quadratic breeding rate is critical in the sense that a BBM 
with breeding rate [3(x) = const • x p explodes in a finite time a.s. if and only 
if p > 2, with explosion in the expected population size when p = 2. When a 
branching inward OU process with quadratic breeding is considered here, a strong 
enough drift with /i > a\/2b can balance the high breeding, whereas any lower drift 
would lead to a dramatically different behavior. o 

Example 11 (Outward OU process with constant breeding rate). Let a 2 ,fi,b > 
and consider 

L := -a 2 A + fix ■ V on R d 



(1 - e - 2 ( Q / CT2 )*) 



exp 
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corresponding to an 'outward' Ornstein-Uhlenbcck process and let (3() = b. As the 
spatial movements have no affect on the branching, the global population grows like 
e* 34 and this is achieved 'naturally' with particles moving freely. This corresponds 
to (L + (3)<j) — b<fi with = 1. On the other hand, the principle eigenvalue is 
\ c = b — (i < b with <f>(x) — const • exp{— {fi/a 2 )x 2 }, it being associated with the 
local, as opposed to global, growth rate. 

After some similar expectation calculations to the inward OU in quadratic poten- 
tial, an upper bound on the process' spread is roughly the same as for an individual 
outward OU particle, that is, we can take a(t) = exp{(l + 5)(fi/a 2 )t} for any 6 > 0. 

Despite the transient nature of the original motion, the spine is an inward OU 
process 

2 V0 _ 1 2 . _ md 



{L+/3-\ c y 



L 



1 , A 

—a A — fix ■ V on 



with equilibrium <j)<t>{x) oc exp{— (fi/o~ 2 )x 2 }. Intuitively, this is the motion that 
maximizes the local growth rate at A c (here its the original motion 'conditioned to 
keep returning to the origin'). We can therefore take ((x) = (1 + e)(a 2 /fi) log a: for 
any e > and hence still find that C( a t) — (l + e)(l + 5)i = 0(t). All the conditions 
required for the strong law to hold are again satisfied. o 

Example 12 (BBM with (3 € C+(R d ) and [3 ^ for d = 1,2). Consider the 
(|A + /3)-branching diffusion where f3 S C+(R d ) and ^ for d = 1,2. Since 
Brownian motion is recurrent in dimension d — 1,2, it follows that A c > and in 
fact, the operator | A + (3 — A c is product-critical and even product-p-critical for all 
p > 1 (see Example 22 in [BT] ). 

We now show how to find a £ that satisfies (ra*) in Definition [U We do it for 
d = 1, the d = 2 case is similar. 

Let b > be so large that supp(/3) c [—6,6] and let M := maxR/3. Recall 
that p(t,x,y) denotes the (ergodic) kernel corresponding to (|A + (3 — A c )^. In 
this example P will denote the corresponding probability. By comparison with the 
constant branching rate case, it is evident that a t := \J2M ■ t is an appropriate 
choice. Therefore we have to find a ( which satisfies that for any fixed ball B, 



lim sup 



p(x,BX(t)) 



1 







together with the condition that ((a t ) — £(V2A/ • t) — 0{t) as n — > oo. 

An easy computation (see again Example 22 in ET]) shows that onl\ [—b, b], 

1 



2A + /3-A, 



\A- S gn(x)-V^c± 



where sgn(x) := x ^ 0. Fix an e and let r±t and To denote the first hitting 

time (of a single Brownian particle) of [—6, b] and of 0, respectively. We first show 
that as t — > oo, 

t(l + ey 



(5) 



sup P x 

b<\x\<t 



r±b > 



0. 



Obviously, it is enough to show that for example 

t(l + e 



Pi 



T > 



STRONG LAW OF LARGE NUMBERS FOR BRANCHING DIFFUSIONS 



11 



where P corresponds to — -\/2A c 4- on [0, oo). Indeed, if W denotes standard 
Brownian motion starting at the origin with probability Q, then 



r > 



*(l + <0 

^/2X 





Yt(i+«) > 


= 




% /2A c 





X 2A, 



t(l 



/2X 



+ W t(i+o > 



W t(i+.) > 



(the last term tends to zero by the SLLN for W). 
We now claim that £(t) := ^^j^ satisfies 



lim sup 



p(x,B,((t)) 



- 1 



= 0. 



(The condition £(a t ) = O(i) is obviously satisfied.) By the ergodicity of p(t,x,y), 
it is sufficient to show that £ satisfies 



lim sup 

*^°°&<|x|<t 



p(x,B,((t)) 



J B #(y) d v 



- i 



= o. 



Let, for example b < x < t. By the strong Markov property at Tb (the hitting time 
of b) and by ©, 



p(x, B 



m _ P (b,B, m -^) 



J B <M>(v) d v 

uniformly in b < x < t. 
Finally, 



lim 

t — >oo 



J B #(y) d v 

p(b,B,0) 



P 



T b < 



i(l 



/2Ac 



/2A7 



Jb ^(v) d v 



= 1 



because pit, x, y) is an ergodic kernel and 

t(l+e) _ 



lim 

£ — >oo 



C(*)- 



/2X 



= lim 



= OO, 



completing the proof of our claim about £. 



Example 13 (Bounded domain). First note that when D is bounded, an important 
subset of V P (D), p > 1 is formed by the operators L+/3 which are uniformly elliptic 
on D with bounded coefficients which are smooth up to the boundary of D and 
with A c > 0. That is, in this case L + (3 — A c is critical (see [P95 , Section 4.7), 
and since <f> and <f> are Dirichlet eigenfunctions (are zero at the boundary of D), it 
is even product-p-critical for all p > 1. Theorem [7] thus applies. 

Although in this case Y is not conservative in D, in fact even Theorem [6] will 
be applicable whenever (in*) can be strengthened to the following uniform (in x) 
convergence on D: 



(0) 



lim sup 

f-> oo xeD , y eB 



p(x,y,C(t)) 



<H>(y) 



o. 



(Note that |AH77j has a similar global uniformity assumption — see the paragraph 
after Q.) Indeed, then the proof of Theorem [B] (which can be found later, in 
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Section |4j can be simplified, because the function a is not actually needed: D an 
can be replaced by D for all n > 1. 

As far as ^ is concerned, it is often relatively easy to check. For example, 
assume that d = 1 (the method can be extended for radially symmetric settings 
too) and so let D = (r, s). Then the drift term of the spine is b + a(log </>)'. Now, 
if this is negative and bounded away from zero at s — e < x < s and positive and 
bounded away from zero at r<x<r + e with some e e (0, s — r), then ([6|) can 
be verified by a method similar to the one in the previous example. The above 
condition on the drift is not hard to check in a concrete example (note that since 
4> satisfies the Dirichlet boundary condition, \og<j) tends to — oo at the boundary). 

If we relax the regularity assumptions on L + [3 then for example <j) is not nec- 
essarily upper bounded, and so we are leaving the family of operators handled in 
|AH77j (see the four paragraphs preceding ((4|); nevertheless our method still works 
as long as L + [3 € V*{D), p > 1 (for the SLLN) or L + (3 G V P {D), p > 1 (for the 
WLLN). ' o 



4. Proofs 

4.1. A spine approach. To establish the L p (P$ x ) convergence of for p > 1 we 
appeal to a, by now, standard techniques that have been introduced to the literature 
by [LPPj and involves a change of measure inducing a 'spine' decomposition. Similar 
applications can be found in [5][BK, EK, HHc] to name but a few. See for example 
[Ev| IE] as well as the discussion in EK] for yet further references. 

It is important to point out that we will need the spine decomposition not only to 
establish the L p -convergence mentioned above but also in the key lemma (Lemma 
ITT]) in the proof of Theorem [5] In both cases, we found the spine method to be 
indispensable and we were not able to replace it by other L p methods. 

Before we can state our spine decomposition, we need to recall some facts con- 
cerning changes of measures for diffusions and Poisson processes. 

Girsanov change of measure. Suppose that Y is adapted to some filtration 
{Gt '■ t > 0}. Under the change of measure 

= g^*) e - K(\c-MY.))ds 

the process (Y, Pj) corresponds to the /i-transformed (h = <j>) generator (L + (3 — 
A c )^ = L + acj)~ 1 V <fi ■ V. Note now in particular that since L + f3 £ V\ (D), it follows 
that (Y, Pj£) is an ergodic diffusion with transition density p(x, y, t) and an invariant 
density </></>. 

Change of measure for Poisson processes. Suppose that given a non- 
negative continuous function g(t),t > 0, the Poisson process (rt,L 9 ) where n = 
{{<Ji : i = 1, ■■■,n t } : t > 0} has instantaneous rate g(t). Further, assume that n is 
adapted to {Q t : t > 0}. Then under the change of measure 

= 2"*expj-^ ff(s)tfej 
the process (n, L 2s ) is also a Poisson process with rate 2g. See Chapter 3 in [JSj . 




dh 2 3 
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Theorem 14 (The spine construction). Let {Tt : t > 0} be the natural filtration 
generated by X. Define the change of measure 



dPs* 



cj){x) <t>{x) 

J~t 



Then, under Ps x , X can be constructed as follows: 

• a single particle, Y = {Y}t>o> referred to as the spine, initially starts at 
x and moves as a diffusion corresponding to the h-transformed operator 
L + a0 _1 V^ • V; 

• the spine undergoes fission into two particles at an accelerated rate 2/3(Y) 
at time t, one of which is selected uniformly at random to continue the 
spine motion Y; 

• the remaining child gives rise to an independent copy of a P-branching 
diffusion started at its space-time point of creation. 



An similar construction for BBM was established in Chauvin and Rouault [CRj . 
See [EK (Theorem 5) or HHa , for example, on how to prove it. 

Remark 15 (The Spine decomposition.). Theorem H4l savs that (X, Ps^) has 
the same law as a process constructed in the following way. A (Y, P^)-diffusion is 
initiated along which (L, 13; D)-branching processes immigrate at space-time points 
{(Y^jCij) : i > 1} where, given Y, n = {{o-j : i = 1, ...,n t } : t > 0} is a Poisson 
process with law L 2 ^( y ). It will often be very useful to think of (X, Pg a ) as being 
constructed in this richer way and it will be convenient to define the natural filtra- 
tion of the spine and the birth process along the spine as Q t '■= cr(Y s ,n s : s < t). 
Note that using the 'spine construction' of (X, P) as discussed in Remark Qj)l we 
can write 

nt 
i=l 

where, conditional on the spine filtration Q t , Wi is an independent copy of the 
martingale Wf started from position Y ai and run for time t — (Tj where cr, is the i th 
fission time along the spine for i = 1, . . . , n t . Remembering that particles off the 
spine behave the same as if under the original measure P and that the martingale 
property gives Ps x (Wf) = <fi(x), we then have the so called 'spine decomposition': 

n t 

(8) E{wt\Qt) e-^V(^) + E e " Ac<Ti ^) 

i=l 

The L p -convergence of the martingale. The a.s. convergence of can be 
complemented with the following result. 

Lemma 16. Assume that L + [3 belongs to V P {D) and that (f3(f> p , 4>) < oo for some 
p e (1, 2]. Then, for x G D, W* is an LP (P$ x )- convergent martingale. Moreover, 
if swp D [3 < oo, then the same conclusion holds if we only assume that p > 1. 

Proof. Pick p so that q = p - 1 G (0, 1], (<p p , 4>) < oo and (f3(j> p , 0) < oo. (If K := 
swp D (3 < oo and we only assume that p > 1, then ((fP, <j>) = {(j) p ~ 1 , <t>4>) < oo implies 
{(jf^) = {(jf- 1 ,^) < oo and {f3<f,4>) < K '{(j?- 1 , #) < oo for all r G (0,p), and 
so we can assume that in fact p £ (1,2].) 
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We adopt an approach similar to the one in |HHa| IHHbj . position Y ai and 
along the spine for i = 1, . . . , rit. original measure P then have the so called Using 
the conditional form of Jensen's inequality, the spine decomposition ([5]) and that 
(u + v) q < u q + v q for u, v > when q £ (0, 1), we find 



{Wff 



(w?y]=E Sx {E[(wfy\g t ]} 



< e Sx {[e (wf\g t )] 9 } 



i=l 



E*[0(y t ) 9 ]+E* 



2P(Y s )4>{Y s y&s 



Call the two expressions on the right hand side the spine term, A(x, t), and the sum 
term, B(x,t), respectively. Since Y has generator L+ a0 _1 (V</>)-V and (4>,<f>) = 1, 
Y is ergodic and Ef(f(Y t )) -> (f4>,<j>). Then 

lim e x " qt A(t, x) = limE^((^(y t ) 9 ) = (cf> p ,4>) < oo 

ffoo tfoo 

for all x £ D. For the sum term note that lim sToo (/3(r s )0(y s )«) = 0) < oo 

and so \\m. t ]oo B(t, x) < oo for all x £ D. By Doob's inequality, is therefore an 
L p -convergent (uniformly integrable) martingale, as required. □ 

4.2. Proof of Theorem [6] along lattice times. The statement that Es x (W^) = 
1 follows from Lemma [T5] The rest of the proof will be based on the following key 
lemma. 

Lemma 17. Fix 5 > and let B CC D. Define 

U t = e- x ° t (<t>\ B ,X t ), 
where c/>\b{x) = <p(x)l^ x ^By Then for any non- decreasing sequence {m n } n >i, 
lim \U/ m +n \ s - E(U( mn+n)s \Fn8)\ = 0, P$ - a.s. 

njoo 

Proof. We will suppress the dependence in n in our notation below and simply 
write m instead m n . Suppose that {Xi : i = 1, iV^s} describes the configuration 
of particles at time nS. Note that we may always write 



(9) 



U, 



(m-\-n)5 



where given J- n g, the collection {{/, 



(i) 



l,...,iV n ,s} are mutually independent 



and equal in distribution to U m s under Pg x . respectively. 

By the Borel-Cantelli lemma, it is sufficient to prove that for x £ D and for all 
e > 0, 

(m-\-n)5 

— EiJJ ^ m+n y\Tn&)\ > e) < oo. 

n>l 

To this end we first note that, 

Ps x (\U(m+n)s - E(U( m+n ) S \T n s)\ > e) < — Es x (\U( m+n y - E(U \ m +n)s\^ : ns)\ P ) ■ 
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Now recall the following very useful result, for example see [B] or |CH| : if p £ (1,2) 
and Xi are independent random variables with E{Xi) — (or they are martingale 
differences), then 



(10) 



E 



<2* ^ElXif. 



Jensen's inequality also implies that for each n > 1, | Yn=i Ui \ P — nP 1 S"=i(l u i| p ) 
and, in particular, \u + v\p < 2P- 1 (\u\ p + \v\ p ). 
Note that 



U, 



s+t 



N t 



i=i 

r(0 _ EVTtWi 



where conditional on jF t , := J7 S — E(Us \Ft) are independent with E(Zi) = 0. 
Thus, by (| 10[) and Jensen, 



i?(|c/ s+t -£;(c/ s+f |^)l p 1^: 



»=i 

i=l 

< 2 p e" pAc *^2 p - 1 £:(|c/ s w | p + £ ; (l^i 4) ri-? r t) -Ft) 



i=l 



< 2 2p e~ pAct ^£(V.. 



{i)\P 



Then, as a consequence of the previous estimate, we have that 

EMI" 



n>l 



(11) 



< 2 2 p £ e - A =" 5 P£; 5:c J2 E sx< KUmsT} 



vi=l 



Recalling the definition of the terms A{x, t) and -B(cc, i) from the proof of Theorem 
ITBT ii) and trivially noting that Ut < Wf , we have 

(m+n) \Fn5)\ P ) 

n>l 

< 2 ^ £ e~^E 5x ( Es Xi 

n>l \l=l / 

< 2 2 ^£ 5:c (^e-^ nS cf>(X i )(A(X i ,md)+B(X i ,mS))\ 

n>l \i=l / 

= 2 2 p <f>(x)e- qXcSn Et (A(Y nS ,m5) + B(Y nS ,m5)) 



n>l 
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where we have used the 'one-particle picture' (equation I|2j0 and the spine change 
of measure at (0. Since the spine Y is Markovian and ergodic under P£, we know 
that 

E* [A(Y n5 ,mS)] = e- x ^ nS Ei(^Y (m+n)s y). 

Denoting m,oo := lini„^oo m n , the latter converges to e _<?Ac " loo ' 5 (0 p , (p) (which will 
be zero if m^ = oo) as n | oo. Recall the assumption that (/3(/) p , <f>) < oo. Similarly 
as before, we have that 

E£ [B{Y nSl mS)} = 2 / e- A ^ s Ef (/3(y s+n5 )0(F s+ „ 5 )«)ds 



which has a finite limit equal to 2 J"™ e _AcS ((3cj) p , </>)ds as rt | oo. These facts are 
enough to conclude that the last sum remains finite to complete the Borel-Cantelli 
argument. □ 

We now complete the proof of Theorem [S] along lattice times. Assume that 
L + (3 &V* for some p > 1. Recall now that 1(B) := J B cj)(y)(j)(y)d(y) < 1 and note 
that, similarly to ([9|), 

E{U t+s \T t ) 

N t N t 
N t N t 

= ]>>- A ^(W(5) + E e ~ Ac( ' +s) / Hy)Q(Xi,y,s)d( y ) 
i=i i=i Jb 

N t . 

= i{B)w? + Y. e ~ Xc(t+s) / 4>(y)Q(Xi,y,s)d(y)=:I(B)wt + e(t,s). 

i=l 75 

Let us replace now i by nS and s by m„<5, where 

m n := ((a nS )/5, 

and £ is the function appearing in the definition of V* . (Although we do not need 
it yet, we note that, according to (iv) in Definitional one nas m n < Kn asn-> oo, 
where K > does not depend on 5.) Then 

E{U [n+mn)S \T nS ) = I(B)wf lS + Q(nS,m n 5) = : I(B)wf lS + T(n). 
Define the event 

A n := {supp(A n<5 ) ^ -D a „J- 
Because of the choice of m n and since 1(B) < 1, we have 

|T(n)| < 

iV„j 

Y,^ n5 cb(Xt)e- x ^ 5 a mn s + \T(n)\l An = e~ A ^ ; a m „ 4 W*j + \T(n)\l An . 

i=l 

Since, according to Definition djiv) , lim n ^oo 1a„ = 0, P— a.s., therefore 
limsup |T(n)| < lim e - A ^ 5 a m „ 5 = P 4 . - a.s., 

n|oo "Too 
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and so 

(12) lim E 5x {U {n+mn)s \T n5 )-{<j)\B,4>dx)Wi =0 P s „-a.s. 

n^oo 

Since Span{0|s, B CC £>} is dense in C+, the result for lattice times follows by 
standard arguments along with Lemma [TTl □ 

4.3. Replacing lattice times with continuous time. The following lemma is 
enough to conclude the convergence in Theorem[6] (see the remark after the lemma) . 
It upgrades convergence along lattice times to the full sequence of times and is based 
on the idea to be found in Lemma 8 of |AH76j . 

Lemma 18. Suppose that for some p > 1, (<j) p , (f>) < oo and for all 6 > it is true 
that for all g S C^(D) and x G D 

lim e-^ nS { 9l X nS ) = (g, 4>)Wt P s=c - a.s., 

n\oo 

then the same result holds when n5 is replaced by t and lim„ by lim^oo . 

Remark 19. Recall that we assumed that C( a t) — 0(t) as t — * oo, and so referring 
to the previous subsection, m n = C( a ns)/S < Kn with some K > which does not 
depend on 8. In fact, by possibly further increasing the function a, we can actually 
take C( a t) = Kt and m n = Kn. Then, from the previous subsection we already 
know that 

hm„ Too e-^ K+1 ^ 5 (g,X [K+1)nS ) = (g,4)W* P S=C - a.s. 
Thus the assumption in Lemma [TSl is indeed satisfied (write 5' :— 8(K + 1)). O 
Proof. First suppose that B CC D and for each x G D and e > 0, define 

B'(x) = {y e B : <j>(y) > (1 + e)-^(x)}. 
Note in particular that x € B e (x) if and only if x 6 B. Next define for each 6 > 

^■ S b( X ) = l{supp(A' t )CB e (a;) for all te[0,<5]} 

and let S,g e (x) = Es x ('E. B x (x)). An important feature of the latter quantity in the 
forthcoming proof is that £g e (x) — > 1b(x) as 5 j 0. With this notation we now 
note the crucial estimate 



e 



-s 



e^ 4 (0| B ,X t ) > 7 f_ T >>-^»V(x i )S£ e (X i ). 



(1 + e)^ 

Note that the sum on the right hand side is of the form ^ where now Ut m + n \§ is 

played by the right hand side above and U^ s is played by the role of (j)(X i )E S ^(X i ). 
Similar L p estimates to those found in Lemma [T71 show us that the estimate flT]) 
is still valid in the setting here and hence 

y~]Es x (\U( m+n )s - E(U( m+n ) S \T n s)\ p ) 

n>l \i=l / 
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However, with q = p — 1, the righthand side can again be upper estimated by 
2 2 P e- x ^E Si (p, X n5 ) = 2 2 " £ e- x ° nS «Et(tf>{Y n5 y) < oo 

n>l n>l 

where the equality follows by equation ([2]), and the fact that the final sum is finite, 
follows by the ergodicity of P!* and the assumption that (<^> p , 4>) < oo. 
We may now appeal to the Borel-Cantelli Lemma to deduce that 



lim 

n^oo 



i=l 



= 



P^-almost surely and hence, using the fact that the Strong Law of Large Numbers 
has been proved already for resequences, 

-\ c t/l\ \r \ c /AC"' 1 



nminfe-^(0| B) X t ) > _ -(<^ e , 
t|oo (1 + ej 

Taking S J. reveals that {4>£, S B ,4>) — > {(j)\ B ,4>) in the lower estimate above, and 
hence subsequently taking e J. gives us 



liminfe-^^lB,^) > 



Recall that this estimate was computed for the case that Bcc fl. Suppose now 
that BCD (not necessarily bounded). Then there exists an increasing sequence 
of compactly embedded domains in B, say {B n : n > 1}, such that U n >i ^ n = 
We may then note that for each n > 1 

liminfe- A ^| s ,X t ) > liminf {<j)\ Bn , X t ) > (^|b„,0)W*, 

tfoo ifoo 

and hence, as n is arbitrary, 



liminfe- Act (0| B ,Xt) > (</>\ B ,<l>)W*, P Sx - a. 



S. 



t]oo 



Now that we have a tight lower estimate for the liminf for arbitrary Borel BCD, 
we shall look at the limsup, also for arbitrary Borel BCD. Using the normalization 
(</>, </>) = 1, one has 

limsupe- A ^| s ,X t ) =Wi -liminf (^\ D \ B , X t ) < {(f>\ B ,$)W+, Ps x - a.s. 

This, together with the liminf result, yields 

lim e- Kt {cp\ B , X t ) = {4>\ B ,j>)Wi, P Sx - a.s. 

f\ oo 

Then, just like for lattice times, a straightforward measure theoretic consider- 
ation shows that <f>\ B can be replaced by an arbitrary test function g G C^(D), 
completing the proof. □ 

4.4. Proof of Theorem E 

Proof. The last part of the theorem is merely a consequence of Lemma 1161 For 
any g £ C+(D) define for each x G D the function h s (x) — Ef(g(£ s )), and note 
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that, uniformly in x and s, the function h s (x) is bounded. Now define Ut[g] 



-A c t/ 



, Xt) and observe that, just as in Theorem [6l one has 



i=l 

where by 

E(U®\g]\r t )=<f>{X i {t))h a (X i (t)). 

Next, note from the Markov property at t and the proof Q of Theorem [6] that for 
fixed s and e > 

Jim Pa. (\U t+s [g] - E{U t+s [g]\T t )\ > e) = 0. 
By the Markov inequality, for each e > 0, 



P Si (^E{U t+s [g\\Ft)-{<l>g,ct>)Wt 



E(U t+s [g]\T t )-(cj>g^)W! 



However, making use of the many-to-one identity, the right hand side above respects 
the following inequalities, 

E S \E{U t+s [g]\Ft)-{<i>gA)W« 



N t 



< 



Hence taking limits as t f oo, and using ergodicity of the spine as well as the uniform 
boundedness of h s (x), we have 

cj>{x) 



-(\h s -(<t>g,<f>)\,<f>cl>). 



Finally, noting that lim^oo h s {x) = ((f>g,(j)} and again the uniform boundedness of 
h s (x) 7 we have by dominated convergence that 



lim lim P,5 

si oo f\oo 



E(U t+s [g]\^t) - {<f>g,4>)W? 



>e < 



e sf oo 



lim \h s - {cf>g, </»)|, #) = 0. 



To finish, we need an epsilon-delta argument. First note that, under the given 
conditions, the martingale converges in the L p norm and hence, in particular, 
converges in probability to its limit. Next, fix e > 0. Then 



Ps x (\u t+ s[g] 



> e 



< P Sx ( \U t+s [g] - E{U t+s [g]\Ft)\ + E{U t+s [g]\^t) ~ {4>g,<t>)W} 

+(6a,6) Wf-Wt 



> e 



< 



Ps x (\U t+s [9]-E(U t+s [g]\F t )\ > 

+P S . C ({<Pg,$)\wf -W+ 



Ps. 
c 



E(U t + s \g]\r t ) - {fafiW* 



> 



> 



^Note that even though Ut is defined differently, we still have martingale differences and the 
key upper estimate of Ut < const ■ wf still holds. 
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Now for each 5 > 0, by choosing s sufficiently large, we can take limits on the left 
hand side above to deduce 



>£]<*. 



So, in fact 



hmP fe (u t+s [g]-{4>9,<t>)Wi 

tjoc V 



lmiP 5x (u t [g]-{ct>g,<j>)Wt >e)=0. 



In particular, taking g = n/cj) for any n e C+(D), yields the convergence (in 
probability) of exp{— X c t}X t in the vague topology of measures. □ 
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